Tidal heating and torquing of a Kerr black hole to next-to-leading order in the tidal 

coupling 
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We calculate the linear vacuum perturbations of a Kerr black hole surrounded by a slowly- varying 
external spacetime to third order in the ratio of the black- hole mass to the radius of curvature of the 
external spacetime. This expansion applies to two relevant physical scenarios: (i) a small Kerr black 
hole immersed in the gravitational field of a much larger external black hole, and (ii) a Kerr black hole 
moving slowly around another external black hole of comparable mass. This small- hole/slow-motion 
approximation allows us to parametrize the perturbation through slowly-varying, time-dependent 
electric and magnetic tidal tensors, which then enables us to separate the Teukolsky equation and 
compute the Newman-Penrose scalar analytically to third order in our expansion parameter. We 
obtain generic expressions for the mass and angular momentum flux through the perturbed black 
hole horizon, as well as the rate of change of the horizon surface area, in terms of certain invariants 
constructed from the electric and magnetic tidal tensors. We conclude by applying these results to 
the two scenarios described above. 
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I. INTRODUCTION 

The Kerr metric describes the external spacetime of a 
spinning black hole (BH) in General Relativity and it is 
the most general, stationary and axisymmetric solution 
to the vacuum Einstein equations [l|-l6| . Real astrophysi- 
cal BHs, however, are not in vacuum, but rather subject 
to vacuum, gravitational perturbations from surround- 
ing objects. These perturbations will result in a change 
of the BH intrinsic parameters, characterized by mass 
and angular momentum flux through the horizon. Since 
this flux could affect the gravitational- waves (GWs) emit- 
ted by BH binary systems through a modification to the 
balance law and, thus, the GW frequency evolution, it 
is important to fully understand the horizon dynamics 
involved. 

The effect of tidal perturbations on GWs emitted dur- 
ing BH inspirals varies depending on the system consid- 
ered [rl-fioj. On the one hand, for quasi-circular, compa- 
rable mass BH inspirals, the effect of these fluxes on the 
GW phase amounts to less than one radian for an event 
in the LIGO frequency band Q. On the other hand, for 
quasi-circular, extreme mass-ratio inspirals, such tidal ef- 
fects can increase the duration of the signal up to twenty 
days in two years for an event in the frequency band of a 
space-borne detector pT] - [T3| . For eccentric and inclined 
EMRIs, this effect could be even larger [l^ - [l6| . 

Poisson 17[ developed a general formalism for calculat- 
ing the linear, vacuum, gravitational perturbations of a 
BH due to a dynamical external universe. In this formal- 
ism, one assumes the external universe is slowly-varying 
and the perturbations it produces can be expanded in 
powers of the ratio of the BH mass to the radius of cur- 
vature of the external universe. This external universe is 
parametrized with the Weyl tensor, which can be written 
in terms of electric and magnetic tidal tensors and their 



derivatives. With these expansions, the equation that 
governs the evolution of linear perturbations, the Teukol- 
sky equation, can be separated and perturbatively solved. 
Its solutions can then be used to compute the fluxes of 
mass and angular momentum across the perturbed hori- 
zon, as well as the increase in the horizon surface area. 

Poisson's work is applicable to both rotating and 
non-rotating BHs through the use of different frame- 
works: the metric formalism and the curvature formal- 
ism. The metric formalism consists of directly perturb- 
ing the Schwarzschild metric and solving the linearized 
Einstein equations. This perturbed metric allows for 
the calculation of the Regge- Wheeler [l^ and Zerilli [l^ 
master functions, which can be then used to calculate 
the horizon fluxes. The curvature formalism consists of 
perturbatively solving the time-domain Teukolsky equa- 
tion (20l - l22j for the Newman-Penrose (NP) scalar ipQ. 
This scalar is then used to obtain generic formulas for 
the horizon fluxes. Although the Schwarzschild case has 
been studied extensively mostly through the metric for- 
malism [l3, [23l - [28j , the Kerr analysis through the curva- 
ture formalism has only been carried out to leading-order 
in inverse powers of the radius of curvature of the exter- 
nal universe [13, [H [H]. 

The purpose of this paper is to solve for the linear vac- 
uum perturbations of the Kerr metric to next-to-leading 
order in inverse powers of the radius of curvature of the 
external universe through the curvature formalism, so as 
to obtain the associated tidal fluxes. We work in the 
curvature formalism because a typical perturbation of 
the Kerr metric, in terms of tensor spherical harmonics, 
does not allow for the separation of the linearized Ein- 
stein equations in all coordinates. Instead, we choose to 
work with curvature perturbations, and more specifically 
with the NP scalar ipo, which will be obtained by solving 
the Tcukolsky equation. 
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We begin by assuming that the external universe 
causes a vacuum perturbation that is slowly-varying. 
The Weyl tensor of the external universe can then be 
decomposed into electric and magnetic tidal tensors that 
depend only on time. After projecting this Weyl tensor 
onto an appropriate tetrad, we obtain an asymptotic TpQ^^ 
that automatically satisfies the Teukolsky equation in the 
asymptotic region. We construct a full ipQ by promoting 
the radial dependence of t^o^^ to arbitrary functions of 
radius. These functions are such that satisfies the 
Teukolsky equation not just asymptotically, but down to 
the BH horizon, while approaching ^q'^^ in the asymp- 
totic region. With ipQ in hand, we can compute the mass 
and angular momentum fluxes, as well as the increase in 
horizon area within the curvature formalism. 

The main results of this paper arc this full NP 
scalar, together with the mass and angular momentum 
fluxes derived from them in terms of certain invariants 
that parametrize the external universe geometry [see 
Eqs. (jlOOp and (jlOip ]. These quantities scale as 



(J) - o{M^/n^) + o{M^/n^) , 



(1) 



(2) 



where M is the mass of the perturbed BH and J the 
magnitude of its spin angular momentum. Similarly, the 
rate of change of the horizon's surface area scales as [see 
Eq. (fT02)l ] 



(A) - oiM^n'^) + o{m^/'rJ') . 



(3) 



Certain assumptions in the small-hole/slow-rotation 
approximation prohibit us from taking the slow-rotation 
limit and comparing our results to Schwarzschild BH 
ones. To achieve this comparison we modify these as- 
sumptions to obtain expressions valid in the slow-rotation 
limit of the curvature formalism. With this, we recalcu- 
late the NP scalar, as well as the horizon fluxes and the 
horizon area increase. We show that the results obtained 
with these new assumptions reduce correctly to the ones 
obtained in the metric formalism when considering lin- 
ear vacuum perturbations of the Schwarzschild metric. 
Therefore, we will here provide two sets of expressions for 
i/'o Enid the tidal fluxes: one set valid for slowly-rotating 
BHs and one set valid for rapidly-rotating BHs. The 
construction of a set of expressions valid uniformly for 
all rotations is left for future work. 

This paper is divided as follows: Section [H] describes 
the curvature formalism and the perturbativc scheme cm- 
ployed; Section IHII computes the asymptotic form of the 
Weyl tensor and the NP scalar; Section ITVl computes the 
exact form of the NP scalar i/'o- Section |V] computes the 
Teukolsky potential ^' and the mass and angular momen- 
tum fluxes; Section |VT] computes the mass and angular 
momentum fluxes in the limit of a slowly-rotating BHs; 
Section IVHI applies our results to certain astrophysically 
motivated scenarios. Section IVHIl concludes and points 
to future research. 



Henceforth, we employ the following conventions. We 
use geometrized units, where G = 1 and c = 1. The 
symbol 0{a) stands for terms of relative order a. Greek 
indices range over spacetime coordinates, while Latin 
indices in between parentheses denote tetrad compo- 
nents. The Einstein summation convention is assumed 
all throughout the paper, where repeated indices arc to 
be summed over unless otherwise specified. 



II. THE CURVATURE FORMALISM AND THE 
SMALL-HOLE/SLOW-MOTION 
APPROXIMATION 

A. The Curvature Formalism 

The curvature formalism is described in detail in [l7j 
and here we only present the basic ideas that are relevant 
to this paper. 

Let us consider a background Kerr BH with mass M 
and spin angular momentum with magnitude J = aM . 
Its perturbations can be characterized by the NP scalar 
'00 J which can be written in terms of the perturbed Weyl 
tensor 5Caf}-yS as 



(4) 



where P and m" are two tetrad 4-vectors. This scalar, 
when computed in the Kinnersley tetrad, diverges at the 
unperturbed horizon. Therefore, it is more convenient to 
express ipo in terms of the Hartle-Hawking tetrad. Using 
ingoing Kerr coordinates {v,r,9,ip), that are related to 
Boyer-Lindquist coordinates (teL, ''bl, 6*61,, 0bl) through 
the transformation 



W = ^BL + 



-dr , 



= <t>BL + a 



A 



-dr , 



(5) 
(6) 



the curvature variable ^ (the so-called Teukolsky poten- 
tial) is then free of divergences on the horizon 17| and 
defined by 



0, V) = -0o(HH) I 



rV'o(K) 



4(r2 + a?Y ^^^^ 
(7) 

Here, '0o(^) and V'o(HH) are the NP scalars in the Kin- 
nersley and Hartle-Hawking tetrad respectively, A 



- 2Mr 



and r± ~ M ± V AP — is the radial 



location of the unperturbed Kerr horizon. The axial sym- 
metry of the Kerr solution allows us to decompose 4* in 
azimuthal modes 



(8) 



that evolve independently. 
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Using the dynamics of the perturbed horizon, we can 
calculate the change in its surface area via [Ol 



.{A) 



r\ + a? 



Using the first law of BH mechanics 



-{A)^{M)-nH{j). 



and the fact that 



(M) 



7n,uj ~ \'-' /m.uj 1 

m 



(9) 



(10) 



(11) 



for any (m, u)) Fourier mode, the mass and angular mo- 
mentum fluxes are [l^l 



(M) = 



4k 



■E 



2k I ( I 1 2) sin 6'd6' 



imnH / (i>™$™ - $™$'") sin^dfl 



(12) 



and 

{j) 



4k 



E(^ 



(13) 



where 



$™(«,0)= r 

^ — cx: 



e-(«-*™"")"'*™(w',6')dw', (14) 
™"«^'*™(w',6')dw', (15) 



are integrated curvatures. In the above equations k = 
(r+ — M)/(r^ + a?) is the surface gravity of the unper- 
turbed Kerr horizon and 51// = i/ (?'+ +a^) is the angular 
velocity of the unperturbed horizon. 



B. The Small-Hole/Slow-Motion Approximation 

For a non-rotating BH, there are only two relevant 
length scales: the mass of the BH M and the radius 
of curvature of the external universe TZ. The assumption 
that the external universe is varying slowly implies that 



M 



(16) 



If we assume that our perturbed BH is in a circular orbit 
around some other object of mass M^^t, then TZ is given 
by the GW wavelength, which leads to 



M 

n 



M 



M + M„t 
where V is the orbital velocity. 



< 1, 



(17) 



The above requirement can be met in two ways: the 
small-hole approximation or the slow-motion approxima- 
tion. In the former, M/M^^t ^ V~^, so our perturbed 
BH is immersed in the gravitational field of a much 
larger BH and the orbital velocity is unrestricted. Of 
course, since V cannot exceed one, this condition trans- 
lates to M/M„^t ^ 1. In the slow- motion approximation, 
V -4: {1 + Mext/M)i/3 so the BH moves slowly and the 
masses of the BHs are unrestricted. Since the masses are 
positive, this condition can be met if one requires that 
y < 1. 

The rotating case is slightly different, since now we 
have two length scales associated with the perturbed BH. 
One of them is again M, while the second is the time scale 
associated with rotation 



1 



+ _ 2M(i + yr^) 



X 



(18) 



where we have defined the dimensionless spin parameter 
X = a/M. The requirement that the external universe 
be slowly-varying implies that both M and T^^t be small 
compared to TZ. The second condition gives 



M 
TZ 



< X, 



(19) 



and since < x < 1 for Kerr black holes, this is a stronger 
constraint than Eq. p6p . This, in turn, changes the defin- 
ing conditions of the slow- motion/small- hole approxima- 
tions into 



M 



M + M„ 



-^'«x- 



(20) 



Thus, in the small-hole approximation, we must require 
that M/Mcxt ^ x/^^; while in the slow- motion approx- 
imation we must have V ^ X^^'^(l + Mo^t/M)^^^ . For 
equal-mass, slowly-rotating BHs in a circular orbit, which 
come into contact with a velocity of roughly V ^ 0.4, 
the above condition puts a limit on how small x can be, 
namely x ^ 0.03. 

The above constraint prohibits us from taking the 
X = limit and comparing our results to those obtained 
within Schwarzschild BH perturbation theory [l7j . The 
slowly-rotating case requires a different analysis, where 
our expansion parameter no longer is forced to satisfy 
M/TZ < X, but rather M/TZ = 0(x). Physically, this 
is because in the x ~^ limit, the rotational timescale 
formally diverges (see Eq. (|18p ). and thus, we cannot 
require that the radius of curvature of the external uni- 
verse be much larger than it. The slow-rotation limit of 
the small- hole/slow- motion approximation is described 
m more detail in Sec. Ell 



III. ASYMPTOTIC FORMS 

Since the perturbations we are dealing with are of vac- 
uum type, in the sense that there is no matter in the 
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vicinity of the background BH, the effect of the exter- 
nal universe is fully described by its Wcyl tensor in the 
neighborhood of the background BH. Our first task is, 
therefore, to determine its asymptotic form through or- 
der 0(7^-3). 

To perform this computation we may ignore the grav- 
itational field of the BH, which at r ^ r+ is negligible 
compared to the tidal effects of the external gravitational 
field. If we imagine that the BH moves on a timelike 
geodesic 7 in the external gravitational field, we may ex- 
press the external Weyl tensor as a Taylor expansion in 
powers of r at fixed v (27j . The Taylor expansion is im- 
plemented formally as 



CafJ-yS — do 



Ca 



(21) 

where Cap-yS is the Weyl tensor evaluated at a point 
X off the world line, while Ca'piy'S' is the Weyl tensor 
evaluated at the advanced point x' on the world line, 
which is linked to a; by a future-directed (from x to x') 
null geodesic. The quantity a{x,x') is half the squared 
geodesic distance between x and x' , a a' is its gradient 
with respect to , and {x,x') is the parallel propa- 
gator from X io x' . 

Following the methods described in Appendix D of [l^] , 
we erect a tetrad (m" , e" ) on the world line 7, which 
we parallel transport off the world line with g°'^, . Ten- 
sors on and off the world line can be decomposed on the 
tetrad, which we indicate for example via CaQbo{x) — 

Cayfjse'^uielu^ , in which = g°'a'<^a and it" = g^a'"^"' ■ 
It may be shown that the frame components of the paral- 
lel propagator differ from the Minkowski metric by terms 
of order and higher, so that they play no role in an 
expansion of the Weyl tensor through order r. Making 
use of ct"' = ru"' - x^e"^' — Eq. (D4) of [13 — we find 
that the frame components of the Weyl tensor are 

Caoboiv, r, 0, 'ip) = Caoboiv, 0) - rCaoboiv, 0) 

+ x'Caoboleiv, 0) + 0(rV7^"), (22a) 

Cabco{v,r,0,'ip) = Caf,co(w, 0) - rCabcoiv,0) 

+ x''Cabc0\e{v,0) + O{ryn^), (22b) 

Cabcd{v, r, 0, = Cabcdiv, 0) - rCabcd{v, 0) 

+ x'Cabcd\e{v, 0) + 0(rV7^^), (22c) 

in which x° ~ [rsin0cos'ip,rsui0 simper cos 9] 
is a Cartesian system defined in the usual way 
from the spherical polar coordinates {r,0,il)), 
CaObo = Ca'^'p'S'-e'e'^ u<' el u^' u"', and Camo\e = 

C'o-'-y'/g'i'ie'e" e\ , with the other projections 

defined in an analogous way. 

The frame components of the Weyl tensor and its 
derivatives evaluated at r = can be expressed in 
terms of the irreducible components provided by the tidal 
quadrupole moments Eah and Bab and tidal octupole mo- 
ments Eabc and Babe which are all symmetric-tracefree 



(STF) tensors in the Cartesian system a;°. The relevant 
relations are listed in Appendix D of , but we repeat 
them here for completeness: 



CaOfco(w,0) 


— ^abi 


(23a) 


Ca6co(w,0) 




(23b) 


Cabcdiv,^) 




(23c) 


Ca060|e(^'7 0) 


^ ^ab\e-! 


(23d) 


CQ6cO|e(^',0) 


— ^abpiJ c|e' 


(23e) 


Cafccci|e(«,0) 


~ ~^abp^cdq^ u 


(23f) 



with the covariant derivatives of the quadrupole tidal 
tensors 

Eab\c = Eabc + ^ {eacpB\ + ^bcpB^a) : (24a) 

4 1. 

Bab\c = -^Babc - g {<^acp£\ + '^bcpS^a) ■ (24b) 

It is a straightforward task to insert the irreducible de- 
compositions within Eq. to obtain an expansion of 
the Wcyl tensor through order r. 
We can now compute the NP scalar 



(25) 



where /" and m" arc members of the null tetrad involved 
in the decomposition of the Weyl tensor. In our compu- 
tation the nuU tetrad at [v, r, 0, V') differs from the null 
tetrad at {v,0,d,ip) by terms of order r^, which can be 
ignored. Therefore, the null vectors are given by = 1, 
1°- ~ n°- ^ x°- /r, and 

m" = — p(cos cos?/)— i sin?/), cosO sinip+i cos?/', — sinfl). 
V 2 

(26) 

With this, we find that 

■00 = -Caobam'^m'' + 2C abcon" m'' ~ Cacbdn^vrfn^m'^ . 

(27) 

An explicit computation of ipo through order r is facil- 
itated by the introduction of the quantities 

"2,0 = i'li + ^22, (28a) 

a2,±i = f 13 T -1^23, (28b) 

a2,±2 = ^11 - ^22 T 2i5i2, (28c) 

p2M = Bii+B22, (28d) 

p2,±i^Bi3TiB23, (28e) 

/32,±2 = Bii - B22 T 2iBi2 (28f) 

and 

Q!3,o = ^113 + £223, (29a) 

a3,±i = fill + ^^^122 T «('S'ii2 + '^^222), (29b) 

a3,±2 = ^113 - ^^^223 T 2zfi23, (29c) 

a3,±3 = fill - 3£'i22 T « (3^112 - ^222), (29d) 
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^/33,0 = ^113 + B223, 



(29e) 



^/?3,±i = Sill + S122 T «(Sii2 + S222), (29f) 



^/33,±2 = 'B113 - ^6223 T 2iSi23, 



(29g) 



^/33,±3 = Sill - 3Si22 T «(3Sii2 - S222), (29h) 

which act as substitutes for the hnearly independent com- 
ponents of the tidal moments. This calculation is also 
aided by the introduction of the spin-weight s = +2 
(scalar) spherical harmonics 



again Eq. (|32)) . We have carried out the calculation both 
ways and found agreement. 

IV. THE NEWMAN-PENROSE SCALAR 

In this section, we compute the full NP scalar ipo asso- 
ciated with a perturbed spinning BH metric. It is defined 
by Eq. (|4]) and in what follows, we use its asymptotic form 
constructed in Sec IIIII to determine the full scalar that 
solves the Teukolsky equation everywhere inside r <C 7?., 
including at r = r_(.. 



A. Teukolsky equation 



2Y^^ -sin6l(cos6lT 1)6^*''', 



and 



2-f3 



(30a) 
(30b) 
(30c) 



(31a) 



2^3^! -sin6'(3cos6'± l)(cos6iT l)e^"^, (31b) 
8 

2y3±2 = i(3cos6l±2)(cos6lT l)^e±2''^, (31c) 



2Y^-^ = -sin6'(cos6'T 1) e 



2 ±3iiA 



(31d) 



These functions are mutually orthogonal, but for our con- 
venience they are not normalized. 

Putting all of this together, we find that the asymptotic 
form of the Weyl scalar is given by 

V'o ~ c = - E 2ir - ^ E (^, v) 

m rn 

+ \rY. 2ir (0, V) + 0{r^/n^) , (32) 



where we have defined 



(33) 



and an ovcrdot indicates differentiation with respect to v. 
The asymptotic NP scalar satisfies the Teukolsky equa- 
tion in the asymptotic region, which we have verified ex- 
plicitly. 

An alternative way to derive ip'^^ goes as follows. 
One starts with the metric of a vacuum spacetime ex- 
pressed in null coordinates (w, r, 6, (j)) in a neighborhood 
of a geodesic world line situated at r = 0; this metric 
is presented in Sec. Ill B of Ref. [l^. One then com- 
putes the Weyl tensor for this metric, and projects it 
onto the tetrad described previously. The end result is 



It is convenient to work with the Fourier transform of 
i/'o, which can be resolved in mode functions that separate 
the variables: 



^0 = E zi^{'^)Rem{r) 2S'"\e)t 



imip 



(34) 



where Rem{r) are free functions of radius, 28^"^ {6) are 
spin- weight +2 spheroidal harmonics and the overhead 
tilde stands for the Fourier transform. The angular func- 
tions reduce to spin-weight s = +2 (scalar) spherical 
harmonics in the limit as becomes time-independent 
(the zero- frequency limit). Notice that all advanced time 
derivatives of tpo simply bring down factors of —iuj when 
working in terms of the Fourier transform. 

This separation of variables allows us to decouple the 
Teukolsky equation in Kerr coordinates into an angular 
equation [H, ^ 



1 



■do (sine 8028 



Aau! cos 



sm 



4m cos ( 



sm 



- 4cot^ + - 4 28'"" = O(cj'), (35) 



and a radial equation [22|, |30| 

A drrR'^'" + 2 [3(r - M) iK] a^i?^™ 



8i{r - M)K 



+ Qiujr — A 



R 



0{lo^). (36) 



Here, E^"^ is a separation constant, we have defined K = 
(r^ + a^) u) ~ am, A = E^"^ — 2amuj — 6, linearized in uj 
and we have set s = +2. These expressions correct a typo 
in Eq. (2.10) of the term explicitly proportional to 
UJ inside the square brackets of Eq. ([55)1 should really be 
2(2s-l)ia;r, instead of -2{2s+l)iujr in Eq. (2.10) of Hi]. 

One can recast the radial equation in a more amenable 
form by transforming to 



(37) 



Doing so, the radial equation becomes 
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x{x + + (6x + 3 + 2i7m) d^W"" + 



4i7m-^^^i4 - + 3)(^ - 2) 



x(l + a;) 



R 



em 



2iio{r+ -r^)\ [x(l + x) + {2p - l){x + p)] + 
£{£ + !) +4 



2{2x + 1) [x{l + :r) + (2p - + p)] 
a;(l + x) 



3(x + p) — imj- 



£(1 + 1) 



R 



em 



(38) 



where we have defined 7 = a/(r_|_ — r_), p = r_|_/(r+ — r_) 
and we have used the eigenvalue 



^ £ (£ ^_ 1) 



8am 
£{£ + !) 



tu + Oiuj^), (39) 



in the w ^ 1 limit for s = +2. Note that this equation 
reduces exactly to Eq. (9.25) in [13] when ui = and 
£ = 2. 

The rest of this section is devoted to finding the solu- 
tions to the angular and radial sector of the Teukolsky 
equation to linear order in lo. 



B. Angular Sector 



rt3±3 v3±3 „t3H' 



(47) 



The £ = 3 modes do not require any 0{uj) corrections. 
These spheroidal harmonics satisfy the angular equation 
in Eq. ([55]) with the eigenvalues of Eq. ([M)) to order 
0{uP'). For future convenience, we will also write the 
spheroidal harmonics as 



where 



2 16 2 

5'45'"9 



1,3m 



. (49) 



Let us first concentrate on the angular sector. Spin- 
weighted spheroidal harmonics satisfy the angular equa- 
tion. In general, they must be solved for numerically, 
but in the limit of small u one can represent them as 
a finite sum of spin- weighted spherical harmonics of the 
same weight [2l| 



^gim^ 2r^"(^»-4aw^ 

i^e' 



2£+l 
2£' + 1 



{£lm0\£'m){£l^20\£'-2) , 

whose eigenvalues were given in Eq. ([39]). The quan- 
tity {jii2mi'm2\J M) arc the usual Clcbsch-Gordan coef- 
ficients. 

In this paper, we are interested in the £ ~2 and £ — Z 
modes only. With our choice of normalization for the 
spin- weighted spherical harmonics, the angular functions 
are given by 



25^° = 


v20 2 3Q 

21 — -auj 21 , 
5 


(41) 


= 


( 2^2^! + gaw 2r-^±i) e^'^ , 


(42) 


^2±2^ 


2Y^^^ - ^aw 2r3±2^ gT2.^ ^ 


(43) 


25^° = 


2^30 


(44) 




^y3±l ^ 


(45) 


53±2^ 




(46) 



C. Radial Sector 

Let us now concentrate on the radial sector. Usually, 
when working to all orders in oj, Eq. (|38p is solved nu- 
merically, after imposing certain boundary conditions at 
infinity and at the horizon. When working to linear or- 
der in uj, however, solutions can be found analytically 
in terms of hypergeometric functions [2l|, [3l|, |32j . When 
expressed in closed form, these solutions are rather com- 
plicated and require the imposition of certain recursion 
relations. For this reason, we will re-derive them here 
and present explicit solutions at 0(w*') and ©(w) for the 
£^2 and £ = i modes separately. 



Expansion of the radial function 



Let us then decompose the radial functions as 



R 



R 



emfi 



(50) 



where R^™--^ and i?^™'^ are assumed to be w-independent. 
With this decomposition, the radial equation to ©(w^) 
becomes 

{x{x + l)dxx + [3 {2x -t- 1) + 2i7m] 
4 



x{l + x) 



[ijm (1 + 2x) 



^{£ + 3){£-2)x{x + l) 



R 



em,o 



= 0. 



(51) 
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Note that this equation reduces exactly to Eq. (9.25) 
in [TtI when £ = 2. Similarly, the radial equation to 
linear order in uj is 

^x{x + l)d^j: + [3 {2x + 1) + 2ijm] 
4 



x{l + x) 
1 



[i^m (1 + 2x) 



{£ + 3)ie ~ 2)x {x + I) 



R 



e.m.1 



T 



£m 



(52) 



to leading order in uj. 

The radial function must also satisfy a regularity con- 
dition at the event horizon. Because R^"^{r) is associated 
with a Weyl scalar "00 constructed with the Kinnersley 
tetrad, which is singular on the event horizon, the radial 
function itself cannot be expected to be nonsingular at 
r = r+. Instead, the relation of Eq. ([7]) implies that it is 
A^i?^™(r) that must be smooth at r = r-|_. In terms of 
the new radial coordinate x, this means that x^i?^'"(a;) 
must be smooth at a; = 0. 



where the source T has nothing to do with a matter 
source, but rather it is constructed from the zeroth-order 
solution and its first time derivative, namely 



T^" = 2i(r+ - r_) I [x{l + x) + {2p - l){x + p)] 
"2(2x + 1) [a;(l + x) + {2p - l){x + p)] 



3(a; + p) — im7 



x{l + x) 
£(f +l) + 4 



R 



Im.O 



(53) 



2. Boundary Conditions 



Asymptotically, we know that ipQ should approach 
Eq. as r 3> r+. Let us then impose this condition to 
see what the boundary conditions on Rgm should be at 
infinity. Fourier transforming Eq. p2p and setting this 
equal to Eq. (f34| we have 



E 

[z2mR: 



(54) 



where we have replaced Z2m ~ ^i^Z2m and we have 
rewritten the spheroidal harmonics as in Eq. (|48p. 

Let us now look at different £ modes separately. 
Clearly, all terms proportional to 2^^'" must be set to 
zero independently from those proportional to 2^"^™ ■ 
Let us then decompose the radial function as in Eq. ([50|) 
with f = 2, so as to separate the term of 0(1) from the 
term of C(w). For the {£,m) = (2,m) mode to 0(1), we 
have 



R' 



1, 



while to 0{ijj) we find 



2mS 



R 



Similarly, for the {£, m) — (3, m) mode we have 



R 



3m, 



(55) 



(56) 



(57) 



3. Radial Solution: m = Case 

Let us first concentrate on the solution when m = 0, 
as this is easier than for general m. To leading order in 
w, we find the solutions 



R 



20,0 



■20,0 



a 



20,0 



i2x + l){6x^ 



6 In 

' 6a; 



1 + x 
1) 



2x2(1 + a;)2 

^30,0 ^ ^30,0 ^ _^ ^30,0 

240a;3 + 130^2 + 120^4 



(58) 



60 (2a; - 
10a; 



l)ln 
1 



(59) 



a;2(l + a;)2 

where Cf'^j'' ^'^d 2'° are constants of integration. 

Before, we argued that the boundary conditions should 
be chosen such that in the asymptotic region r » r+ , ipo 
behaves as in Eq. p2l) . In particular, this implies that 
the boundary conditions should be given by Eqs. (|55p 
and ([57)) . Asymptotically expanding the two solutions, 
about a; = 00, we find 



^20,0 

^30,0 



Ci"'" + 0(x-5), 
2Cf '"a; + 0(1) . 

Imposing Eqs. ([55]) and ([57)) . 

1 . ^ 

2^ 



Cf 



1 , 



CI 



30,0 



-(r+ -r_) 



(60) 
(61) 



(62) 



The selection of C|° and Cf^ is based on the require- 
ment that x^R^"^{x) must be smooth at a; = 0. Expand- 
ing the solutions about a; = 0, we find that 



R 



20,0 



' — ^ -f - -I- 6 In X 
2a;^ 



R 



30,0 



C 



30,0 



1 



X 

12 

X 



0{l) 
60 In a; + 0(1) 



(63) 
(64) 



In order to eliminate the logarithmin terms, we must 
choose C2°'" = = C^'^ . To leading order in w, the 
solutions then become 



^20.0 ^ 1 



R 



30,0 



a f 1 



(65) 
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Let us now move on to the radial solution to linear 
order in cj for ^ = 2. The 0{uj) piece of the ^ = 3 solution 
is not needed as it would contribute to 0(7?,""') in ipQ. 
The £ = 2 solution contains a homogeneous piece plus 
an inhomogeneous piece. The homogeneous piece is the 
same as that of Eq. ([58| . Including the inhomogeneous 
solution, we find 



R 



20. 1_ 



61n| 



1 + x 



(2x + l)(6a;2 + 6x- 1) 



2x2(1 +x)2 



+ ln(x) 1^ h4i7a + 2iM j - ln(.T + 1) 

+ —X - iM(\ + X r - ^ 1 + a^r 

07 D 



4z7a 



(66) 



where constants of integration. 

Let us now impose the boundary conditions. Asymp- 
totically expanding the £ = 2 solution about x = 00, we 
find 



i?20,i ^ ^ 2iA/ln(x) + Cf + 0[ln(x)] 

37 



(67) 



The solution automatically satisfies the boundary condi- 



tions in Eq. ([57)1 . We are free to choose ' as we wish, 
and one possible choice is to eliminate the constant term 
as one asymptotically approaches infinity via 



/-.20,1 * 



(68) 



Asymptotically expanding the £ ~ 2 solution about x 
0, we find 

^20.1 0^20,1 / 
i?204 _ „^ + £^+in(x) ( 6Cf ' V2*Af + - + 4*7a 



2x2 



7 



(69) 

The above solution possesses two types of divergences at 
the horizon; a logarithmic and a polynomial one. The 
requirement that defined in Eq. be smooth on 
the horizon forces us to eliminate the logarithmic one by 
setting 



20.1 



2170 ia 
~i 67 



iM 



(70) 



^. Radial Solution: m 7^ Case 

Let us now return to the general to 7^ case. The 
solutions at leading-order in w, ie. to Eq. (|5T|). are 



R 



2m,0 



A2rnX~'Hl + x)-2f(-4, 1; -1 + 2iTO7; -x) , 



(71) 



R 



3m, 



yl3„x-2(i + a;)-2i^(_5, 2; -1 -I- 2iTO7; -x) , 



(72) 



where F(a, b; c; z) is the hypergeometric function and we 
have defined 



^2m = ~ g"^7 (1 + *"^7) (l + 4m272) , 



(73) 



= -T^rna (1 + imj) (l + Am^Y) (3 + 2imj) 
ioU 



(74) 



These solutions have been normalized such that R^ 



1 and R 



3m, 



r as X 



oo and such that they have a 



well-defined Taylor expansion about x = 0. Note that no 
natural logarithm terms appear in this solution, although 
clearly they both diverge as x~2 as x — >■ 0. 

At linear order in uj we again focus on the i = 2 mode, 
as the £ = 3 mode would contribute at (D{TZ~^) to ipQ. As 
before, the solution will contain a homogeneous piece plus 
an inhomogeneous one. The homogeneous piece will be 
identical to Eq. (|7ip . and thus, it is redundant and should 
be discarded via the appropriate choice of integration 
constants. 

One is then left with the inhomogeneous solution, 
which can be obtained with the method of variation of 
constants. To implement this we need both homogeneous 
solutions of the radial equation. The first is given by 
Eq. ([7T|) . while the second is 



Ro (x) = 1 + - 

' X 



2i7m 



(75) 



and the inhomogeneous solution is given by 



R 



2771, 1 



-R^ (x) 



R. 



,2m, 



(x) 



Rt 



2m,0 



(07^'" (0 



di, (76) 



where W is the Wronskian associated with R^™' and 

r,2m,0 

The homogeneous solution R^"^'^ deserves further dis- 
cussion. When solving for the leading-order in to radial 
function we did not present this solution because it fails 
to produce a x'^R^™'{x) that is smooth at x = 0. We 
therefore set the integration constants such that this so- 
lution would not contribute. To linear order in w, how- 
ever, one needs to use both R^'^ and F^^'^ to construct 
the inhomogeneous solution. 

The first integrand in Eq. (|76p is a polynomial in ^ 
and can be integrated immediately. The second term 
leads to a sum of hypergeometric functions, which can 
be manipulated with the identity 

a{h-c)zF{a + l,b;c+l;z)=c{c-l)F{a- l,fe;c- l;z) 
-c[c-l + {a-b)z]F{a,h]c;z), (77) 

which can be derived from the Gauss relations for hy- 
pergeometric functions. By repeated application of this 
identity, the second term can be expressed in a form that 
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3 2 

features 2-^1 [1, 1; 3 + 2im7; —a^] only. The inhomogeneous ^ ^ l_ ^ ^ ^ ^ ^ ^ x 

solution then becomes 87 (a; + 1)^ * '^^ (x+l)^ (.t + 1)^ 



,2™,i_ x2F(l,l;3 + 2i7m;-x) 



i? = r where we have defined 

(i — ^mj [x + 1) 



^ _ (1 + K„)(l + 2k [IOk^ - (5 + 18p)Km + 3p(8p - 1)] 

^ _ 1284, + (164 - 264p)Kf„ + (4 - 552p + 576p^)Kf„-(41 + 294p - 720p^)4^-(9 - 30p - 72p^)K„,-36p(2p - 1) 

^ 36ip(l - 2Km)K,n 

(79b) 

^ 16^ + (20 - 24p)4, - 24p(3 - Ap)^^ - (5 + 42p - 96p')kI - (1 ^ Wp + Sp')^^ - 4p(2p - 1) ^^^^^ 

2ip(l + 2k„i)(1 - 2Km)K„, 
^ 32^^ + 324„ + (4 - 120p + 192p^)4 - (8 + 72p - 144p^)^^ - 3(1 - 2p)(l - 4p)«^ - 6p(2p - 1) 

6ip(l + Km)(l + 2Km)(l - 2k„)k„ 



in which k,„ = 1771,7 ^.nd we have set the coeflicients 
multiplying the homogeneous solutions to zero. Since 
the hypcrgcomctric function behaves as ln(a;)/a; when 
X ^ 1, in the limit a; — >■ 00 the above solution behaves 
as i?^™'^ and satisfies the boundary condition in 

Eg. ([57]) automatically,. 



D. Summary of Results 

The Fourier transform of the full NP scalar can be 
written as 

jn 

-uY, {R2m,l 2^2" + a b'"' R2mfi 2^3™) 

rn 

- ^3mi?3m,0 2^^™ , (80) 



where the radial functions are given by Eqs. ([7T|) . ([7^ . 
([75)) and the 6^"' coefficients are given by Eq. (jlHl), with 
the (un-normalized) spin-weighted spherical harmonics of 
Eqs. O-dSID- 

Note that although the full asymptotic behavior in the 
asymptotic region TZ ^ r ^ r+ is 



r2m 



■00 ~ - X! -^2171 
m 

m ^ 



(81) 



the third term is subdominant relative to the fourth term. 
Therefore, the asymptotic behavior of Eq. ((80|) is identi- 
cal to Eq. ((51)) . 



V. THE TEUKOLSKY POTENTIAL AND THE 
FLUXES 

In this section, we compute the Teukolsky potential as 
defined by Eq. ([7]). In terms of the quantity x, defined 
in Eq. (|57)) . the Fourier-transformed potential ^ is given 

by 



*(w, 



x^{l + xf 



A[{x + pY+^ 



212 



V'o(K) 



(82) 



Evaluating the right hand side at a; = and making the 
decomposition 



(83) 



we find for the m 7^ case, 

77i7r^ (i — 77?7)(1+ 47^771^) 



«'™(cj,0) =-- 



2880p4Af2 



[3052™ 2y'"(e,0) 



(3 + 277777)z3m 2>^'"(e,0)] 



77772r3 (i-7777)(l+4727n2)52m52" 2>"^"(6',0) 



96p5M2 

and for the 777 = case, 

*520 2Y^\e, 0)(47^+72 + 2pM + r+)rl 



(84) 



^^{LO,e) =LU- 



192A/2p5 



(85) 
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In order to return to the time-domain and express ^ as a 
function of v, we must compute the inverse Fourier trans- 
form of the above equation. One can see that 9) is 
given by the right-hand side of Eqs. and (|55|) . with 
the replacement Z2m 'iz2m/uj. 

The integrals in Eqs. ([14]) and (fTS]) can be evaluated 
perturbatively to the appropriate order. Assuming that 
\]/™ varies slowly on a time scale r = (\i'™/vl/'")~i ^ Jl 
(since the time dependence of 4* comes from external 
universe perturbations) we can write for the m ^ case, 



K — imfl 



1- 



1 



H 



{mVLH)n 



-O 



O 



1 



(86) 



When m ~0, Eq. 
can see from Eqs 



(87) 

87|) clearly cannot be used, but as we 
([T^ and (1131), we will only need 
in this case, and for this we can use Eq. (|86|. Therefore, 
for our purposes, the above equations are valid for all m. 

Let us discuss the expressions for $5:' further. The 
function outside the square brackets has terms of 
0{M/'R?) (proportional to Z2m) and 0{]VP/n^) (pro- 
portional to z^jn and -22?™) i while the second term in the 
square brackets is oi 0{M/TZ). Therefore, combining this 
with the leading-order terms that scale as 0{M/'R?), 
one finds terms of 0{M^ /Ti?). The next-order terms 
scale as 0{AI'^ /TZ^) and the lowest-order corrections that 
these induce will be of 0{M^ /TZ^), which can thus be ne- 
glected. 

Substituting the expression for 4'™ in Eqs. (|86)) and 
(|571). we find 



{i — 7to)(1 -f 47^m^) 



96M2(-K 



Z27, 



ze'™^^«7r^(i - 7m)(l -f 47^771^) 



(89) 



and 



$0 ^ e-"^(4r^7^+2pM + r^)ri 

plus terms that contribute at 0{M'^ /TZ'^). 

The canceling of subleading terms in the above equa- 
tions must be done carefully and keeping in mind that 
<i>™ in the end has to be squared and averaged over. Ini- 
tially, $™ has terms that go as 



The terms in the first line come from the ones in 
the second line come from 'I'™ and those in the third 
line come from 5'™. When squaring the above quantity, 
the two lowest-order terms will come from the squar- 
ing of the Z2m 2^^™ term and the cross term between 
Z2m 2l"^™ and ia™ 2>"^". They will be of 0{AP /W^) 
and 0{M^ /TZ^) respectively. All other terms can be ne- 
glected to the order of approximation considered here. 

Substituting Eqs. (I88l)-(l90l) into Eqs. ([HI) -([131), ex- 
pressing K and flu in terms of M and a, and using the 
dimensionless parameter we find 



2xM 



45 2^ ™ ' 

m/0 



m 



and 



45 



8Af6(cr + 1) 

45" 



2, r (^2 + 3m^x^) 



m 



X (5(22mi2m)) - 2xO-(3?(22m-22m)) } , 



(93) 



where a = \/l — ^2. Notice that this agrees exactly with 
Eq. (9.32) in Finally, using Eq. ^ or the first 

law of BH mechanics, we can calculate the change of the 
horizon area 



E 



45 yr 

128.x^Af<'(. + l)^^^,^^,^,^ 



45 

m Z2m)) 



m/0 



— {'^{z2mZ2m)) 

a 

(94) 



Some comments about the above results are here in 
order. First, one might notice that the above results do 
not depend on the £ = 3 part of ^f™. This is because 
we have here squared $ to obtain the fiuxes. Any ^ = 3 
terms in result in quantities that scale as 0{M'' /TZ^) 
in ( j) and as 0{M^/TZ^) in (M), and they are beyond 
the order considered here. Second, one might also notice 
that the above results do not depend on the second term 
of Eq. (Uni). This term is of 0{M/n^) and it is propor- 
tional to 2^^™- In principle, after the squaring it could 
couple to the leading-order term in Eq. (|80p and give 
another 0{M^ /TZ^) contribution. However, this term is 
proportional to i?2in,i and as a; —> it behaves as 1/x. 
From the above asymptotic behavior, it is obvious that 
this term has no contribution to ^' ~ x'^'iI)q as a: — > 0. 

We can rewrite these expressions using Eqs. ([28l) . ([29]) . 
P3p and the invariants 



•Pi ^ Z2m 


\^2rn 
2 J 


f (23m 2^''" + Z2m 2^'") 


El -- 


_ c cab 
^ ^ab^ 5 


Bi 


= BabB'^' , 


(95) 


+ Z2m 


^y2m ^ 


- {Z3m 2^3™ + Z2„ 2^3") 


E2 = 




B2 


= Babs'B^^s' , 


(96) 


+ Z2m 


2^2™ H 


- O (M^/TZ^) . (91) 


E3 ~ 


{Eabs'^S^'f , 


B3 




(97) 
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E5 = epgc£Pj',s''s''s', B5 



(98) 
(99) 



where s° = (0, 0, 1) is the direction of the BH spin. Doing 
so we find 



(M) 



2M^X 
45 



[-4 (3x' + 1) {Ei + S4) + Ihx^E^ + B5 



(100) 



45 45 

x{-8ax(i;i+Bi)+3ax(x'+4)(£;2 + S2) -3ax'(i;3 + S3)-4(l + llx')(-B4 + S4)+3x'(13 + 2x')(£;5+B5>} ■ 

(101) 



and 



(^) - ^^^[8(1 + 3x^)(i^i + B,) - 3(4 + 17x^)(^;2 + B,) + 15x^(i?3 + B,)] + ^^^^^'^^ + 



45cr 



45 



|8x(^i + Bi) - 3x(x' + 4)(^2 + B2) + SxHEs + B3) + \ [32x'(S4 + S4) - 6x'(4 + x^)(E^ + ^5 



(102) 



VI. SLOW ROTATION 

Equation ([T^ clearly implies that the x = limit is 
inaccessible to our previous analysis. This means that 
Eqs. (jlOOp and (|10ip in the x = limit will not reduce 
to Eqs. (8.38) and (8.39) in This is more clear if 

one considers the assumptions involved in the evaluation 
of Eqs. ([14]) and ([T5|) . As discussed in Sec. |Vl we can 
approximate those integrals with Eqs. ([55]) and ([57)1 by 
assuming that the exponentials vary rapidly compared to 
^f™, which is justified if M /Ti. ^ x- Then, each integra- 
tion by parts brings out a term that is 0{TZ^^) smaller 
than the previous one. This method, however, breaks 
down when the exponentials vary slowly compared to 
^P™, as in the ~ X ^ hmit, which forces Eq. ((87)) 
to diverge. 

The above discussion suggests that the small x limit 
requires a special treatment. In this limit, the timescale 
associated with 51h — x/(4A/) is comparable to TZ: 
X = 0{AI/TZ). The horizon fluxes can be computed by 
expanding the Teukolsky equation in powers of x and 
treating Mw = 0{M/TZ) as a quantity of first-order in 
X- In this section, we calculate specifically the rate of 
change of the BH's surface area, which can be given by a 
simple expression; expressions for (M) and (J) are more 
complicated, and must be written in terms of integrals of 
the tidal tensors. 

We follow the general methods of Sec. IIV[ but use x 
instead of uj as an expansion parameter. Because au = 
O(x^), it follows immediately from Eq. that 

2S'"'{e) = 2Y'"'{0,O)+O{x^), (103) 

with angular eigenvalues given by i? = £{£ -f 1) + 0(x^)- 



r 



The Teukolsky equation then reduces to 

(r^ - 2Mr)^ + Q(r - M)^ - ix2M 
dr^ dr 

r-M / Ar2 \ 3 Ar 



Ar2 

r- — m 

4A//2 

R = 0{x' 



dR 
dr 



(104) 



in which A is defined by w = Re- introducing x ~ 

(r - r+)/{r+ - r_) = r/{2M) -1 + O(x^) and following 
the same procedure as in Sec. IIVI and fVl we find 



1 

12 



M(imr!H^2rn + ^2™) 2i"2"(^, «) + O(x') , 



(105) 

in the time domain. 

The calculation of (A) requires that we first obtain 
('y, 0) via Eq. (fTi)) . Wc evaluate the integral by taking 
advantage of the fact that both e"""^!'' ^nd *'"(«', 6*) 
vary over a time scale TZ that is much longer than ^ 
AM . Integration by parts returns 



— e 



(106) 



up to terms that are smaller by a factor of 0{M/TZ) ~ 
0(x). Inserting this within Eq. (jH]), discarding correc- 
tions of order x^ and defining the new invariants 



Eq — Sab^"" , Bq — BabB"" 



(107) 



gives 



87r 



(A) 



16 

45' 



AnliEi + Bi) - 6nUE2 + B2) 



12 



-4nn{Ei + Bi) + {Ee + Be)\. (108) 
In the limit x — > 0, this becomes 

^{A) = ^M'{Ee + Be), (109) 

which agrees with the Schwarzschild expression in [13] ■ 
When JIh — x/(4-A'/) is large compared with TZ~^, we 
obtain instead 



VII. APPLICATIONS: BH IN A CIRCULAR 
BINARY SYSTEM 

In this section, we apply the main results of Sec. |V] 
to the case where the perturbed BH is in a circular bi- 
nary with another BH with mass M^^^. We first consider 
the case where the orbital velocity is small (slow-motion 
approximation). We then tackle the case where the per- 
turbed BH is small (small- hole approximation). 



2(^1 3(^2 + ^2) , (110) 



A. Slow-Motion approximation 



which agrees with Eq. (|102p in the regime x ^ 1- 

When the tidal field is in a state of rigid rotation with 
angular frequency 57, the Weyl scalar tpo depends on ip 
and V through the combination ip — flv only. This implies 
that the tidal moments Sab and Bab satisfy the relations 



(111) 



in which are constant amplitudes. It is then easy to 
verify that with tidal moments of this form, the invariants 
satisfy the identities 



Eg = 2n^{2Ei - 3E2 



(112a) 
(112b) 



Similar identities are satisfied by the magnetic tidal ten- 
sor invariants. With these, it follows that Eq. (|108l) re- 
duces to 



Stt^ ' 45 ^ 



fin) \2{Ei + Bi) - 3{E2 + B2) 

(113) 

when the tidal fields are rigidly rotating. 

For rigid rotation, (A) gives rise to separate expressions 
for (M) and (J). We rely on the first law of BH mechan- 
ics, {k/8tt){A) = (M) — fl-ii{J), and the rigid-rotation 
relation (Af) = il{J), to obtain 



32 r 

(M) = —M^n(n - nil) 2{Ei + Bi) - :i{E2 + B2) 

45 L 



and 

{j)^-M'{n-nu) 



(114) 



2{Ei + Bi)-'i{E2 + B2) . (115) 



These expressions agree with the results in Sec. IX F of 
171 ] in the regime X'^^- According to our conventions, 
the sign of Vt is measured relative to the BH rotation. 
Assuming that the tidal invariant within square brack- 
ets is positive, we have that the BH gains mass when 
> and > r^n; otherwise the BH loses mass (the 
tidal dynamics is superradiant). We also have that the 
BH gains angular momentum when f2 > f^n, and loses 
angular momentum when f2 < fin. In all circumstances, 
the tidal dynamics produce an increase of surface area. 



Let us first assume a slow-motion approximation, 
where V ^1. If the BH is in a circular orbit with radius 
6, the orbital angular velocity and the relative orbital 
velocity are given by [25| 



n = e 
V -- 



Mt 



Mt 



(116) 
(117) 



respectively, where e = L- s = ±l,?7 = AIM^^JM!^ is 
the symmetric mass ratio and Mt = M + M^^^ the total 
mass. Then, the non-vanishing frame components of the 
tidal tensors are given by [25| 



-A£ii+£22) 



At, 



263 



M 
2Mt 



2(^11-^22)—^^ 



(118) 

cos 2f2t, 
(119) 



£12 



Bi3 



23 



3M^ 
263 

63 
3M^ 
63 



1 



A/,,, + 3A/t 



2A/t 

■V cosVLt + 0{V^) , 
■Vsinm + 0{V^) , 



v^ + o{v^) 



sin 2Q,t , 
(120) 
(121) 

(122) 



to the appropriate order. 

Substituting Eqs. (fTT8)) -([T22 |) into Eqs. pOO| and (fTOTj) 
we find 



^ . , , Af 

^''^--r^[M-T 



x(l + 3x')V' 



15 



6 2( M . 



17 



M 



(23x^-H6)-(l + 3x^)^-(l + 3x^)(^)' 



16 

5 ' \Mt 



(a + l)(l + llx')Fi' + 0(l/^^), 

(123) 



13 



and 



M 



2' ' Mt 



5 ' \Mt) 



(12 + 51x')-(4 + 12x2) 



M 
Mr 



xv 



14 



16 



M 
Mt 



Mt{(t + 1)(1 + llx')'l^'^ + OiV^"^) . 

(124) 



The third term in Eq. (|123p is of higher order that we 
are allowed to keep here. For circular orbits, however, 
(M) = 0,{J), and, thus, we can obtain this term from 
the third term in Eq. p24p . The results agree with [l3| 
to 0{TZ~^). Using the first law of BH mechanics, we find 
the change in the horizon area 

5 \AIt ) <y 



167r 



M 



i3 \Mt, 
10247r 



Mt 



(12 + 51x')-(4+12x') 



M 
Mt 



v2 



-677 



E \ Mt^^^^V'' + 0{V'^). 
Mt J <7 

(125) 



We can repeat the same calculation in the slow- 
rotation limit, as derived in Sec. IVII We find 



(M) 



f.,^(^) MT{n-n,)v^^. 



/ h 32 2 ^ M 



and 



{A) 



10247r 



-e?7 



M 
Mr 



M. 



2 (f2 - f2if)^^i5 



(126) 
(127) 

(128) 



where = x/4M and Vt is given by Eq. (|116p . The 
above equations are valid only up to the order allowed 
by our calculations, given by Eqs. (|123p - (|125p . Notice 
that in the slow-rotation limit, the tidal fluxes are smaller 
by a factor of 0{V^) than when the background BH is 
rotating fast. 



B. Small-Hole approximation 

If the BH is in a strong external gravitational fleld, we 
must let V be unrestricted, and, instead, assume that 
M ^ Afoxt. Assuming the external spacetime is that 
produced by a Kerr BH, and the perturbed background 
Kerr BH moves on a circular orbit of radius b in the 



equatorial plane, the angular and orbital velocities are 
given by 



sgn(^) 



and V = 



(129) 



respectively. In the above equations, S is the dimension- 
less spin parameter of the external BH, with its spin and 
mass given by Jext = SM^^^. 

Since now the external gravitational field is strong, we 
have to calculate the tidal tensors directly from Eqs. ([25)) 
and p4p . The non vanishing components of Sab and Bab 
have been calculated in [2^ and they are given by 



1 



V 1 - 4S1/3 + SE^V^ 



£2 



2A-f^ 



1,„ „ , 3^6 1-2^2^22^4 
2 - ) = - 2Mi: 1 - 3^2 + 2Sy3 2$ , 



£12 = 

Sl3 = 



3^6 1_ 21/2-^22^4 ^ 

2^2,1-3^2 + 25^3 ^^"2$, 



3sgn(S)V^ (1 - Sy)Vl - 2^2 _^ 52^4 



Ml^ 



1 - 3^2 + 2Sy3 
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where 



3sgn(S)V7 (1 - Sy)Vl- 2^2^^21/4 



M2, 



1 - 3V^2 _^ 2Sy3 



(130) 

(131) 
(132) 

sin$ , 
(133) 

cos<i> , 
(134) 



^ ^ sgn(^)y3 

A4,t 



(135) 



Substituting them in Eqs. pOOp and (jlOll) we find 



4 , , / Af 



(136) 



and 



M 
5 



:Sgn(S)((7 + l) 



A/ 



(137) 



where 



1 - 2^2 + S2^4 



[(8 + 9x')s2 V' 



(1 - 3^2 + 2Sy3)2 

+ (6x^ -^)SV^~{A + 27x^) +A + , 

(138) 



1 — 2T/2 I "2t/4 

r. = } ^rJ.... [(6x^ ~ 49x^ - 8) -2 



(1 - 3V^2 + 2Sy3)2 



14 



- {I2x^ - lOx^ - 8) EV^ + {6x^ + 83x^+4) - 44%^ - 4] 

(139) 

The above results agree with ^ to 0(7^"2) ^he 
leading-order term in Eqs. (|136p and (|137p for S = 
(external universe described by a Schwarzschild BH) also 
agrees with Eqs. (9.52) and (9.53) in [l3| and with Eqs. 
(4a) and (4b) in [2g| . The change in the horizon's surface 
is 



16n ( M 



5 \M,, 



I 2 



M 



32^ 



(J 

' X('T + 1) ,5 



F^-^Ta, (140) 



where = Fi + 

As explained in the previous subsection, when we 
work in the slow-motion approximation, the leading- 
order term of Eq. (fTSe]) is calculated from Eq. ([T00| . The 
second term in general is not within the reach of our ap- 
proximations, but since, again, the binary system is in 
rigid rotation, (Af) = SI/ (J), where flf is the angular fre- 
quency defined in Eq. (jl35p . Then, in the slow-rotation 
limit of Sec. IVIl we find 



(M) = ^sgn(S) (^^y MU^-nnWTo, (141) 



/ h S ( M 



l/^'To, (142) 



and 



{A) 



2567r 



sgn(^) 



n 



(143) 

where is given by Eq. (jl29p . Fq = F2(x = 0) and, as 
before, the order up to which the above equations are 
vahd is the same as in Eqs. (fT36)) . (|T37)) and (|T40)) . As 
before, observe that in the slow-rotation limit the tidal 
fluxes are suppressed by ©(Af/Mcxt) relative to the mod- 
erately fast rotating case. 



VIII. CONCLUSIONS 

We studied how the properties of a rotating BH are 
modified when it is subject to external vacuum pertur- 
bations, extending previous calculations to order TZ~^ in 
the slow- motion/small- hole approximation. We mainly 
used the NP scalar ipo to describe these perturbations, 
where '0o must satisfy the Teukolsky equation. 

The asymptotic behavior of ipo depends on the exter- 
nal universe and it is constructed as follows. Working 
under the assumption that the external universe induces 
a small perturbation on the background spacetime and 
it varies slowly (that is, we expand all quantities in pow- 
ers of M/TZ, where M is the mass of the BH and TZ 
is the radius of curvature of the external geometry), we 



parametrized the Weyl tensor of the external spacetime 
through time-dependent symmetric and trace-free tidal 
tensors £ab, Bab, £abc and Babe- Then, the asymptotic 
form of "00 can be obtained by Taylor expanding the Weyl 
tensor to the appropriate order in TZ^^ and expressed in 
terms of the above mentioned tidal tensors. Then, pro- 
jecting the Weyl tensor onto a tetrad, one obtains the 
asymptotic form of ipo . 

This asymptotic NP scalar suggests a decomposition 
of the Fourier-transformed "00 into multipole modes. The 
decomposition introduces some unknown radial functions 
that can be determined by solving the Teukolsky equa- 
tion, thus allowing tpo to also be valid close to the un- 
perturbed horizon, outside the asymptotic region. The 
time dependence of this full scalar is still given by the 
tidal tensors, while its angular dependence is given by 
spin- weighted spheroidal harmonics. 

In order to obtain the fluxes across the horizon, we cal- 
culate tpo on the (unperturbed) horizon, but this quantity 
diverges in this limit if ipo is computed in the Kinnersley 
tetrad. We therefore compute i/jq in the Hartle-Hawking 
tetrad, which is well-behaved at the horizon, and use it to 
compute the mass and angular momentum fluxes across 
the horizon, as well as the change in horizon area. 

No external universe structure is specified when ob- 
taining these results. The only restriction imposed on 
the external universe is that it be slowly-varying and 
its tidal effects on the background BH be small. Once 
the external universe is specified, for example through a 
post-Newtonian metric, the tidal tensors could also be 
determined by asymptotic matching [33l - l35j . 

To conclude, we applied our results to an astrophysi- 
cally motivated scenario: a binary system composed of 
two rotating BHs. We assumed that the external universe 
perturbations are due to a Kerr BH and calculated the 
fluxes of mass and angular momentum and the change 
of horizon surface are of the background BH for two dif- 
ferent cases: a small BH immersed in a large external 
universe and a slow-motion orbit. 

The calculation of the mass and angular momentum 
fluxes across the BH horizon is a definite step towards 
understanding the perturbed dynamics around rotating 
BHs. These fluxes are necessary to self-consistently cal- 
culate the evolution of the orbit, and thus, the emitted 
gravitational radiation. However, in order to fully de- 
scribe the perturbed spacetime, we need its metric, which 
could be a possible avenue for future research. The per- 
turbed metric close to the event horizon of the perturbed 
BH can be constructed from the NP scalar through the 
Chrzanowski procedure [1^, |36| - |39| and it would allow us 
to study the geometry of the perturbed horizon [23, HI] . 

By asymptotically matching this Kerr perturbed met- 
ric to a spinning binary PN metric that is valid far from 
the background BH, we can also obtain a full metric 
that describes the dynamics of the entire spacetime (ssl - 
[ssl . l40j . Such a metric could provide initial data for nu- 
merical relativity simulations [U \^ , or it could serve 
background metric for accretion disk studies [43| . 



15 



ACKNOWLEDGMENTS 

We thank Luis Lchncr for helpful discussions while 
working on this paper. NY acknowledges sup- 



port from NSF grant PHY-1 114374 and NASA grant 
NNX11AI49G, under 00001944. EP acknowledges sup- 
port from the Natural Sciences and Engineering Research 
Council of Canada. 



[1] 
[2] 
[3] 
[4] 

[5] 
[6] 
[7] 

[8] 

[9] 

[10] 

[11] 
[12] 
[13] 

[14] 
[15] 

[16] 

[17] 

[18] 

[19] 
[20] 
[21] 

[22] 
[23] 



R. P. Kerr, |Phys.Rev.Lett.| ll, 237 (1963) 

S. Hawking, LGen.Rel.Gr avj 1, 393 (1971) 

S. Hawking, Commun.Math.Phys.l 25, 152 (1972) 

W. Israel, Phys. Rev. " 164, 1776 (Dec 1967), 

http : //link^aps . org/doi/10 . 1103/Ph ysRev . 164 . 1776, 

W. Israel, Commun. Math.Phys. 8, 245 (1968)" 

B. Carter, IPhys.Rev.Lett., 26, 331 (1971) 

E. Poisson and M. Sasaki, [Phys.Re v. D51, 5753 (1995), 

|arXiv:gr-qc/941 2027 [gr-qc]l 

K. ^n, 'Phys.Rev. D64, 104020 (2001), 
arXiv:gr-qc/0107080 [gr-qc] 

H. Tagoshi, S. Man o, and E. Takasugi,|Prog Theor.Phys.| 

98, 829 (1997), |ar Xiv:gr-qc/9711072 [gr-qc] | 

H. Fang and G. Lovelace, |Phys. Rev. Dp f27 124016 (Dec 

2005) 

S. A. Hughes, Phys.Rev.' D64, 064004 



A. Hughes, Phys.Rev 
,arXiv:gr-qc/0104041 [gr-qc] 
K. Martel, Phys. Rev. D69, 044025 



(2001), 
(2004), 

|gr-qc/0311017| 

R. H. P rice and J. T. Whelan,|Phys.Rev.LettT| 87, 231101 
(2001), 'arXiv:gr-qc/0107029 [gr-qc]" 

A. Pound and E. Poisson, Phys.Rev. D77, 044013 (2008), 
arXiv:0708.3033 [gr-qc] 

N. Yunes, A. Buonanno, S. A. Hughes, M. Cole- 
man Miller, and Y. Pa n, |Phys.Rev.Lett.| 104, 091102 
(2010), |arXiv:0909.4 263 [gr-qc]' 

Buonanno, S. A. Hughes, 
M. C. Miller, and 



N. Yunes, A. 
Y. Pan, E. Barausse, 
W. Throwe, Phys. Rev. D 83, 044044 (Feb 2011), 
|http : //link . ap s . org/d oi7l0 ■ 1103/PhysRevD . 83 . 044044 1 

E. Poisson, Phys.Re^ DTO^ 084044 (2004), 
|arXiv:gr-qc/0407050 [gr^ 

T. Regge and J. A. Wheeler, [Plys. Rev. 108, 1063 (Nov. 
1957) 

F. J. Zerilh, Phys. Rev. D 2, 2141 (1970) 

S. A. Teukolsky, Astrophys. J. 185 , 635 (1973) 

W. H. Press and S. A. Teukolsky, Astrophys. J. 1 185, 649 

(1973) ' 

S. Teukolsky an d W. Press,|Astrophys.J.| 193, 443 (1974) 
E. Poisson, Phys.Rev.Lett. 94, 161103 (2005), 
arXiv:gr-qc/0501032 [gr-qc]. 



[24; 
[25; 
[26; 

[27 

[28; 
[29; 

[30 
[31 

[32; 

[33 

[34; 
[35; 

[36 
[37 

[38' 
[39 

[4o; 

[41 
[42; 
[43; 



K. Martel and E. Poisson, Phys. Rev. D71, 104003 



(2005), |gr^/0502028| 

S. Taylor and E. Poisson, [Phys.Rev.| D78, 084016 (2008), 
arXiv:0806.3052 [gr-qc]' 
S. Comeau and E. Poisson, 



Phys.Rev.] D80, 087501 



(2009), |jarXiv:090 8.4518 [ gr-qc] 
E. Poisson and I. Vlasov , ,Phys.Rev.| D81, 024029 (2010) 
|arXiv:0910.4311 [gr-qc]| 



Poisson, and R. Massey, Cl ass . Quant . GravT] 

D73, 024010 



I. Vega, E 

28, 175006 (2011), arXiv:1106.0510 
N. Yunes and J. A. Gonzalez, Phys. Rev 
(2006), ,gr-qc/0510076. 

S. A. Teukolsky, |Phys' Rev. Lett:] 29, 1114 (Oct. 1972) 
S. Mano, H. Suz uki, and E. Takasugi, |Prog.Theor.Phys.| 
95, 1079 (1996), | arXrv:gr-qc/9603020 [gr-qc] | 
Y. Mino, M. Sasaki, and T. Tanaka, 
Pr og.Theor.PIiy s.Suppl. 128, 373 (1997), 
|arXiv:gr-qc/9712 056 [gr-qc] 
N. Yunes, W. Tich y, B. 
mann(2005), gr-qc/05 03011 

N. Yunes and W. Tichy, |Phys. Rev.| D74, 064013 (2006), 
arXiv:gr-qc/0601046 

N. K. Johnson-McDanieh N. Yunes, W. Tichy, 

and B. J. Owen, Phys.li^ D80, 124039 (2009), 
|arXiv:0907.0891 [gr-qc] | 



J. Owen, and B. Brueg- 



L. Chrzano wski, |Phys. Rev.| Dll, 
Ori, IPhys.Rev.l D67, 



2042 (1975) 
124010 (2003), 



[Phys.Rev. [ 

'arXiv:gr-qc/0207045 [gr-qc]"" 
R. M. Wald, Ph ys.R ev.Lett. 41, 203 (1978) 
L. S. Kegeles and J. M.' Cohen, [Phys. Rev. | D19, 1641 
(1979) 

L. Gallouin, H. NaJsano, N. Yunes, and M. Campan- 
elh(2012), 'arXiv:1208.6489 [gr-qc]' 

T. Chu, Numerical simulations of black-hole spacetimes, 
Ph.D. thesis, California Institute of Technology (2012) 
G. Reifenbe rggr and W. Tichy, .Phys.Rev. D86, 064003 
(2012), arXi v:r205.5502 [gr-qc] | 

S. C. Noble, B. C. Mun dim, H. Naka no, J. H. Kro- 
lik, M. Campan elh, et aL, [Astrophys.J^ 755, 51 (2012), 



[arXiv: 1204. 1073 [astrcHph.HE]| 



